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ABSTRACT 

Nebular spectra of Supernovae (SNe) offer an unimpeded view of the inner region of 
the ejecta, where most nucleosynthesis takes place. Optical spectra cover most, but 
not all of the emitting elements, and therefore offer only a partial view of the products 
of the explosion. Simultaneous optical-infrared spectra, on the other hand, contain 
emission lines of all important elements, from C and O through to the Intermediate 
Mass Elements (IME) Mg, Si, S, Ca, and to Fe and Ni. In particular. Si and S are 
best seen in the IR. The availability of IR data makes it possible to explore in greater 
detail the results of the explosion. SN 2007gr is the first Type Ic SN for which such 
data are available. Modelling the spectra with a NLTE code reveals that the inner 
ejecta contain ~ IM© of material within a velocity of w 4500 km s~^. The same mass 
of ^^Ni derived from the light curve peak (O.OTGM©) was used to power the spectrum, 
yielding consistent results. Oxygen is the dominant element, contributing O.SM©. 
The C/0 ratio is < 0.2. IME account for ~ O.IM©. This confirms that SN2007gr was 
the explosion of a low-mass CO core, probably the result of a star of main-sequence 
mass « 15M0. The ratios of the Ca ii lines, and of those of Fe ii, are sensitive to the 
assumed degree of clumping. In particular, the optical lines of [Fe ii] become stronger, 
relative to the IR lines, for higher degrees of clumping. 

Key words: Supernovae: general - Supernovae: individual: SN 2007gr - Radiation 
mechanisms: thermal 



1 INTRODUCTION in the ejecta, for bo th SNe la and core-collapse SNe (e.g. 

u c 1 1 ■ 1 ■ u J J- iMazzaH et al]|2008al lbl. respectively), 
bupernovae (SNe) are the fanal explosion marking the end ot 

the life cycle of different types of stars. Massive stars (M > properties of the prog enitor st ar may be inferred 

SMq) explode when their core can no longer synthesise new through direct imaging (e.g. [wt^ lioogl ) or searches 

elements and undergoes collapse. Thermonuclear SNe (SNe for the surviving co mpanion in the case of SNe la (e.g. 

la) are the complete explosive destruction of White Dwarfs |Ruiz-LaDuente"^[2004h . However, the only method that 

reaching the Chandrasekhar mass limit, most likely via a^- ^^^-^^ jj^^^^jy ^1^^ gj^ ^-^^^^ ^^^^ ^^^^ transparent 

cretion in a binary system ( |Hillebrandt fc Niemeyeij | 2000D ^ remnant is nebular-epoch spectroscopy. Soon 



although other chanriels have been proposed ( |Pakmor et al.| ^^^^^ explosion, the SN nebula is still optically thick, and 



[2010|;| Finketal. | |201(]| ). There are a number of ways to mea- ^^^^ ^^^^^ 1^^^^^ observed. About 6 months to 

sure the properties of SNe and infer those of the progenitor ^^^^ ^^^^ explosion, however, expansion makes the 

stars: classically, the analysis of SN light curves yields m- ^^^^^^ sufficiently thin that it becomes essentially transpar- 
formation about the mass ejected in the explosion and the radiation. In this phase the gas is heated by collisions 

explosion kinetic energy iArnctt .1982). Combined studies ^-^^^ f^^^ particles produced in the thermalizatfon pro- 
of the light curve and the spectra allow a more precise es- ^^^^ ^1^^ positrons produced in the decay of 
timate of these parameters, as well as of the abundances 56,^^ sep^^ ^^^1^ U^^g emission. As optical depth 

effects are negligible, the energy that is deposited is imme- 
diately re- radiated, whic h makes it easier to estimate the 
* E-mail: mazzali@mpa-garching.mpg.de ^^Ni mass (|Axelrodlll980l ). This phase can be used to study 
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the inner layers of SNe, yielding information on the d etails 
of the explosion, both for SNe la l|Mazzali et al.|[2007al ) and 
core-collapse SNe, especially those that have lost the oute r 
layers (Type Ib/c, iMaeda et ahlbOOSl : iMaurer et al.ll20ld ). 

One major advantage of the nebular phase is that the 
luminosity emitted by a certain ion depends, albeit indi- 
rectly, on the mass of the ion itself. A reliable estimate of 
the emitting mass requires then two basic elements: obser- 
vational coverage of the largest possible number of lines of 
different elements and a reliable calculation of the ioniza- 
tion and excitation state of the gas. For the latter, NLTE 
calculations are usually performed. The former requirement 
implies essentially observing as wide a spectral range as pos- 
sible. Traditionally, because of the faintness of SNe in the 
nebular phase and the easier availability of optical detectors, 
in almost all cases only optical spectra of SNe in the nebular 
phase were obtained. However, several important line emis- 
sion features are located at IR wavelengths. This includes 
elements such as Si, which is very important for SNe and 
has no strong emission lines in the optical, so that its mass 
cannot be derived from optical spectra. Also, the ratios of 
different lines of the same element can be use d to determine 
the d ensity of the emitting gas (e.g. Ca ii, iLi fc McCravl 
1 19931 ). (Fe ii. lLeloudas et alj|2009l ). 

Only a very small number of nebular-phase IR spec- 
tra of SNe is available, and these are typically SNe l a 
spectra (e.g. ISpvromilio et al.ll2004l : iMotohara et al.ll2006h . 
Among core-collapse SN, there are late-time I R spectra for 
the pecuhar SN IIP 1987A (IFassia et al.l 120021') . a nd for the 
highly reddened SNIIP 2002hh l|Pozzo et al.ll2006l ). The first 
stripped-envelope SN for which late-ti me IR spectra hav e 
been obtained is the Type Ic SN 2007gr l|Hunter et al.ll2009^ . 

SNe Ic show neither H nor He in their spectra 
(|Filippenkolll997l). and come with a very diverse range of 
obse rvational (iMatheson et al.l |200G| ) and physical proper- 
ties (|Mazzali et al.l 20061 ) . These depend on the properties 
of the progenitor star and the details of the explosion, rang- 
ing from ejected masses of ~ IM©, ^®Ni masses of Q.IMq 
and kinetic energ ies of ^ 10^^ erg for "normal" events (e.g. 
ISauer et al to ejected masses of ~ IOM0, ^^^Ni masses 

of ~ 0.5Mq and kinetic energies of ~ 5 ■ l O'^^ erg for GRB- 
connected hypernovae (Mazzal i et al.ll2006l ). 

SN2007gr had normal luminosity and spectral appear- 
ance, but was characterised by the unusual presence of 
carbon lines in the early-tim e spectra, indicatin g that C 
was not fully burned t o O ([Va lenti ct al. 2008). A sim- 
ple light curve analysis l|Hunter et al.. .2009 ) suggests that 
SN 2007gr ejected ~ 2Mq of material with kinetic energy 
Ek ~ 1 — 4 10^^ erg, which would make it a rather ordi- 
nary SN Ic. High-resolution radio observations of SN 2007gr 
revealed emis sion from an expandin g source (jParagi et al.l 
I2OIOI . but see ISoderberg eTalT l|2010l )). This expansion was 
interpreted as evidence for material moving at mildly rela- 
tivistic velocities in what seems to be a bipolar jet. This is 
at odds with the existing picture where all SNe associated 
with emission from relativistic material also show large ex- 
pansion velocities in the early phase {v ~ 0.1c), which links 
the relativistic material to the bulk of the SN ejecta, and 
are characterised by large values of the ejected mass and 
the explosion energy (Mazzali et al. ,2006. ). 

We use the nebular spectra of SN 2007gr, including one 
simultaneous optical-IR spectrum obtained more than one 



year after the explosion, to determine the properties of the 
inner ejecta. 

In the following, in Section 2 we briefly review the data 
we use, in Section 3 we discuss the nebular model used for 
the calculations. Section 4 presents the results of one-zone 
models. Section 5 deals with the models obtained using a 
density profile, and in Section 6 our results are discussed. 



2 THE DATA 

Nebular-phase spectra of SN2007gr were obt ained at vari- 
ous epochs, as outlined in lHunter et al.l l|2009 l). Earlier opti- 
cal spectra have epochs of 103 to 158 days after maximum. 
These spectra are not fully nebular. The same is true for var- 
ious IR spectra obtained at similar epochs. At day 375 after 
maximum a simultaneous optical-IR spectrum was obtained. 
Since the optical and the IR spectra were obtained within 
one day of one another, spectral evolution and relative cali- 
bration are not an issue, unlike data at earlier epochs. The 
very advanced phase of this later spectrum ensures that it 
is fully nebular. The sp ectrum was present ed in Figures 7 
(optical) and 12 (IR) of lHunter e\~sl\ (|2009l ). 

Here we use the very late optical-IR spectrum at 376.5 
days after maximum for our calculatiouii^, and verify results 
using the optical spectrum at 158 days after maximum. 

SN 2007gr exploded in a crowded region in the spi- 
ral galaxy NGC 1058. In particular, SN2007gr is flanked 
on either side by a relatively bright stellar association (see 
iHunter et al.ll2009|. Fig. 1) and se veral fainter sources within 
150 PC (see lCrockett et al.ll2008l . Fig. 1). At the late epochs 
under consideration, when the SN is faint in comparison with 
these adjacent objects, particular attention must be paid 
to the estimation of the background flux. In order to esti- 
mate the contribution of these sources we extracted spectra 
obtained in the vicinity of SN 2007gr. The resulting back- 
ground contribution clearly depends on the choice of loca- 
tion at which the flux was estimated. We experimented with 
several extractions along the spatial profile and chose the 
background contribution that best matched the underlying 
slope in the SN continuum. We then scaled this background 
to the apparent SN continuum and subtracted it from the 
SN spectrum (no strong continuum is expected from the 
SN at late epochs). The c ontaminating sources (within 150 
pc) are probably clusters (|Crockett et al.| [ 200 8:) with spec- 
tral energy distributions that peak at blue wavelengths, so 
while some contamination affects the entire spectrum, it is 
expected to be particularly severe at bluer wavelengths. 



3 THE MODEL 

We modelled the spectra of SN2007gr using a code which 
computes the energy deposition from the radioactive de- 
cay of ^^Ni and ^''Co in the SN ejecta, uses the deposited 
energy to estimate ionization and excitation, and balances 
this with gas cooling via line emission in a NLTE scheme 



^ This spect r um w as incorrectly labelled as +375.5 days in 
iHunter et al.l ||2009|) . who did not take into account that 2008 
was a leap year. 
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l|Ruiz-Lapuente Sz Lucy 19921 ). The method foUows that de- 
veloped bv lAxehodI (|l980l ). 

The code aUows the use of clumping, expressed in terms 
of a clumping fac tor Ci which is defined as the inverse o f 
the filling factor (|Li fc McCravlll993l : iMaurer et ahlboid ). 
This assumes that the gas is distributed uniformly, but re- 
sides in clumps. The effect of clumps is to increase the lo- 
cal density, and hence to favour recombination. Clumping 
also reduces the population of the excited levels, which may 
be important in the context of the r atio o f optical to IR 
fiux, as sugge s ted b y iLeloudas et aD (|2009h . For SNelb/c, 
iMazzali et al.l l|200lh found that a typical clumping = 10 
yields best results, while for SNela clumping is not re- 
garded to be necessary if a detaile d model of the explosion 
is used (e.g. iMazzali et al.l l2008al ) . The need for clumping 
in SNelb/c may be the result of an inaccurate treatment of 
ionization. 

Impr ovements in the treatment of ionization 
l|Maurer fc M azzaii| I2OI0I ) make the code more reli- 
able, in particular for SNe where Ni and Fe do not dominate 
the ejecta. These improvements were introduced in the 
context of including a treatment of oxygen recombination. 
A simultaneous improvement of the treatment of ionization 
was also necessary in order to make the calculation of the 
recombination rate more reliable. Forbidden line emission is 
infiuenced by this new treatment only weakly. Results with 
the old and new codes are shown and compared in Section 
4. 

In a one-zone version, the code assumes that the gas is 
distributed with uniform density and abundances within a 
sphere of radius R = v x t, where v is the outer velocity 
of the expanding nebula, as determined from the width of 
the emission lines, and t is the time elapsed from explosion. 
This takes advantage of the fact that in SN ejecta any given 
parcel of matter is characterised by a velocity which remains 
constant as long as the ejecta are not slowed down by in- 
teraction with circumstellar material. Results obtained with 
this model are presented in Section 4. 

In a more sophisticated version the SN envelope is di- 
vided into uniform radial shells. Deposition of 7-rays and 
positrons is computed using a Montecarlo scheme, and emis- 
sion is computed within each shell. This version can be used 
if a detailed model of the explosion is available, which de- 
scribes the run of density with ejecta velocity, or if one is 
sought for from fitti ng the profile of the nebular lines (e.g. 
iMazzah et "al]|2007bl 'l. Results obtained with this model are 
presented in Section 5. 

For the distance and reddening to SN2007gr we use 
the same values (d=9.3Mpc, i. e. distance rnodulu s p — 
29.84 mag; E(B-V)= 0.092) as in lHunter et~al] (|2009l ). 



4 ONE-ZONE MODELS 

For the nebular spectra of SN 2007gr we used first the 1-zone 
code to define the general properties of the ejecta. Here we 
show 1-zone model results obtained first without and then 
with the oxygen recombination module and compare them. 
We then applied the shell model to improve the results and 
the fits. These are discussed in the next section. 
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Figure 1. SN2007gr: one-zone models for the spectrum at day 
172 (29 Jan 2008, black) with the old (blue) and the new version 
of the nebular code, which includes an improved treatment of 
ionization and oxygen recombination (red). 



4.1 One-zone models, old code 

4.1.1 Day 172 

The first nebular spectrum we model was obtained on 2008 
Jan 29. We assign it an epoch of 172 days after explosion, 
using an arbitrary risetime of 13.5 days, which is typical 
for SNeIc and is consistent with the risetime deduced by 
iHunter et al.|[2b09i ') . 11.5 ± 2.7 days. 

The strongest lines in the spectrum are typical of SNe Ic: 
[O i] 6300, 6363 A, Mg i] 4570 A, Ca 11] 7291, 7324 A, the 
Ca II IR triplet near 8600 A, and Na I D near 5900 A. Ca 11 
H&K may be seen at the blue edge of the spectrum. The 
spectrum can be fitted reasonably satisfactorily assuming an 
outer velocity of 5000 km s~^ and — 10, which is typical for 
SNeIc (Figure 1, blue line). The ^^Ni mass wa s kept fixed at 
0.076 Mq , as indicated by the light curve peak (iHunter et al.l 
I2OO9I ). This value yields fairly good results. The mass en- 
closed within 5000 km s~^ is ~ I.4M0, oxygen being the 
main constituent ('^ IMq). 

While most emission lines are reproduced, and identi- 
fied in Figure 1, the main shortcomings of the model are 
as follows. The line near 4050A is not reproduced. This is 
probably [S 11] 4069, 4076 A, which is predicted to be present 
but is much too weak when compared to the data. The low 
ionization degree of sulphur (~ 5% in the models with the 
old treatment of ionization) is probably responsible for the 
weakness of the emission line. The broad base of [O l] 6300, 
6363 A is not reproduced. This is the effect of adopting a 
sharp outer boundary to the nebula. Shell models will im- 
prove this (see Section 5). The emission line at ~ 7700 A IS 
not reproduced. This is most likely the recombination line 
O I 7773 A, which is addressed in the next sub-section. An 
alternative possibility is [S l] 7723 A. This is indeed predicted 
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by the model as the strongest hne of sulphur, but is much 
too weak. The second strongest sulphur line in the model is 
[S ll] 4068, 4076, which is too weak for a S mass of O.OSM©. 
The third strongest feature is [S l] 4507, 4589 A, which is 
completely swamped by the strong Mg l] 4570 A line. In or- 
der to reproduce both the line near 7750 A and that near 
4050A as due to sulphur, the S mass must increase to an 
unreasonable O.65M0, confirming that the line near 7700 A 
is O I 7773. Additionally, cooling by S emission weakens the 
other lines, so that all other masses must increase, to a total 
of ~ 2.4M0. Such a larg e S mass fraction seems to contradict 
nucleosynthesis results (jWooslev fc Weaveilll995l ). 

The line at 8700 A is too weak. This is a blend of Ca ll 
IR triplet and [C l] 8727 A. Increasing the strength of the 
[C l] line is not an option, as the wavelength of the line 
does not match that of the emission peak (assuming spher- 
ical symmetry). Also, the strength of the other carbon line, 
[C l] 9825, 9850 A, is already overestimated. Therefore, a 
weak Ca ll line is the likely cause of the discrepancy. This 
may indicate a diff erent degree of clum ping than adopted in 
the calculation (see lLi fc McCravIl 19931 ). We however refrain 
from investigating clumping here, since this may be better 
done using the optical/IR spectrum of day 390, which has a 
broader wavelength coverage, and the new code, which in- 
cludes oxygen recombination and an updated treatment of 
ionization. Finally, the line at 9300 A, which is mostly due 
to [Co ll] 9340, 9345, is not correctly reproduced. This may 
be caused by insufficient background subtraction. 
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Figure 2. SN2007gr: one-zone models for the spectrum at day 
390 (3 Sept 2008, black) with the old (blue) and the new version 
of the nebular code, which includes an improved treatment of 
ionization and oxygen recombination (red). 



4.1.2 Day 390 

The later spectrum is fully nebular, and it is therefore to 
be expected that our code can reproduce it better than 
the earlier one. Additionally, it covers both the optical and 
the IR essentially simultaneously, making it possible to nar- 
row down the range of acceptable parameters. Many of the 
strongest features are the same as on day 172. One of the 
key features for modelling is the [Fe ll] line near 5200 A. This 
is composed mostly of lines at 5159, 5262, 5273, and 5333 A. 
Other strong lines include the structured blend near 8700 A, 
which is due to a mix of a weak Ca ll IR triplet, [C l] 8727 A, 
and [Fe ll] lines, the strongest being at 8617 and 8892 A (Fig- 
ure 2). The [Fe ll] 8617 A line is actually the strongest of 
all these emissions, which explains why the observed profile 
peaks between Ca llIR and the [C l] line. Uncertainties in 
the strength of the [Fe ll] lines may be responsible for the 
poor reproduction of the corresponding emission at day 172. 

Moving to the IR, several other features are visible 
which allow relative comparisons. They are labelled in Fig- 
ure 3. One of the most useful is [Si l] 1.604, 1.646 fi, which is 
the strongest silicon line, blended with [Fe ll] 1.599, 1.644, 
1.664 jU. Lines of [Fe ll] are also seen near 1.25 n (strongest 
contributors to the feature being lines at 1.249, 1.257, 1.270, 
1.279, 1.321 fi), but the spectrum is very noisy there. Other 
lines of IME are also present in the IR. 

Since the flux in the IR lines of iron is basically de- 
termined from the flux in the optical lines (despite some 
influence of clumping) the abundance of Si can be deter- 
mined from the line near 1.65 /i. This line has contributions 
from both [Si l] 1 608, 1.646 /i and [Fe ll] 1.644, 1.664 
iHunter et al.|[2009l 'l. Given the abundances adopted in our 
model, Fe and Si contribute roughly equally to the observed 
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Figure 3. SN2007gr: one-zone models for the spectrum at day 
390 (3 Sept 2008, black) with the old (blue) and the new version 
of the nebular code, which includes an improved treatment of 
ionization and oxygen recombination (red). 
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Table 1. Results of One-zone Models. 



code V. epoch vel Mass Te log(ne) M(C) M(0) M(Ne) M(Mg) M(Si) M(S) M(Ca) M(56Ni) 
[days] km s'^ Mq K cm-^ Mq Mq Mq Mq Mq Mq Mq Mq 



old 


172 


5000 


1.40 


5100 


7.4 


0.25 


0.99 


6.0e-5 


1.250-3 


0.006 


0.023 


0.053 


0.076 


new 


172 


5000 


1.83 


4870 


7.7 


0.25 


1.40 


6.0e-5 


l.OOe-3 


0.006 


0.070 


0.019 


0.076 


old 


390 


4000 


1.04 


3560 


6.4 


0.05 


0.81 


2.7e-4 


2.80e-2 


0.006 


0.023 


0.048 


0.076 


new 


390 


4000 


1.28 


3430 


6.6 


0.14 


0.99 


3.5e-4 


3.25e-2 


0.009 


0.024 


0.013 


0.076 



emission. Fitting any difference between the observations 
and the synthetic [Fe ll] feature as due to [Si ll] is actually 
the only way to derive the abundance of Si from the nebular 
spectra. The abundance of sulphur can be determined from 
the strong [S l] line at 1.1 /i. This is composed of lines at 1.08 
and 1.13^, which are much stronger than the sulphur lines 
which are seen in the optical. Unfortunately the region of the 
spectrum near 1.1 is very noisy. The ratio of optical to IR 
[Fe ll] emission lines is affected by the degree of clumping in 
the ejecta. However, since we have an independent estimate 
of the ^®Ni mass from the light curve, we can use that value 
(O.O76M0), which basically fixes the amount of iron at day 
390, and modify the filling factor to obtain a best match to 
the optical-IR Fe emission. Any additional flux in the 1.6 /i 
line can then be attributed to silicon, while at 1.1 ^ sulphur 
can be assumed to contribute. 

Using the one-zone model, the spectrum can be best 
matched using an outer velocity of 4000 km s~^. This is sig- 
nificantly less than at day 172. Clearly the outer regions of 
the nebula are becoming tr ansparent, and the emissio n lines 
become narrower (see also iTaubenberger et al.ll2009l ). This 
evolution indicates that the density structure in the ejecta 
is important in shaping the nebular emission lines, and it 
suggests that the approximation of a one- zone model is not 
ideal. As in the model for day 172, we used a ^^Ni mass of 
O.O76M0. A clumping factor = 10 yields the best ratio of 
optical to IR flux. At this epoch the model contains a mass 
of « IM0, and it is again dominated by oxygen (~ O.SMq). 
Other elements have smaller mass. In particular, the Si mass 
required to reproduce the 1.6/^ line is only O.OOGMq. This is 
a very small mass, especially considering that other IME are 
more abundant (e.g. S, Ca, Mg). Still, this model reproduces 
most features. The main shortcoming may be in the [Fe 11] 
line near 5200 A, which may suggest that a somewhat larger 
*^Ni mass may be required. However, background subtrac- 
tion may be uncertain here, since the shape of the feature is 
reproduced reasonably well, as is the ratio of other Fe lines 
in the optical and IR. Interestingly, this ratio increases for 
increasing clumping factor. The higher density in the clumps 
favours recombination and a greater role of coUisional pro- 
cesses, which support the population of the excited levels. 
Otherwise, at the low temperatures that characterise the SN 
nebula at very advanced epoc hs, IR lines sh ould dominate 
(the so-called IR catastrophe. lAxelrodlll980h . We present a 
model with more clumping in Section 5. 

Although the mass of elements such as O and Ca is 
consistent at the two epochs, others change significantly. In 
particular, the Mg mass in the later spectrum is much larger 
than in the earlier one and the same is true, to a lesser ex- 
tent, for Na. Clearly these elements are not very well con- 
strained by the model, probably because of an incomplete 



treatment of ionization. Quite likely, both Mg and Na are 
highly ionized at early epochs and become less ionized at 
later phases. Given the low ionization potentials of Na I 
and Mg I, the high ionization in the early phase is proba- 
bly caused by photoionization and possibly thermal electron 
coUisional ionization, neither of which are treated by the 
code. An estimate of the masses of Na and Mg based on the 
neutral species at times when ionization is still significant 
may therefore yield incorrectly large values if this ionization 
is not taken into account. Mass estimates obtained around 
390 days should be more reliable than at 172 days. 

Conversely, the C mass is much smaller in the later spec- 
trum. Although this estimate is uncertain also because of the 
poor fit of the 8600 A emission, to which [C l] contributes, 
the reduction that is observed may indicate that C is prefer- 
entially located at velocities between 4000 and 5000 km s~^. 
This is not unreasonable given the structure of a stellar 
core: helium seems to "bottom out" at ^ 8000km s^^ in 
SN2007gr, while C is o bserved at v ~ 500 - 8000 km s"^ 
also in the early phase (|Hunter et al.|[2009l . Fig. 14). 



4.2 One-zone models, updated code 

We have recently introduced a number of updates and im- 
provcmen ts to the origi nal code. These are described by 
Maurcr & Mazzalil (|2010l V In particular, we have introduced 
recombination of O I and revisited the treatment of ioniza- 
tion. The former improvement was necessary in order to re- 
produce the O I 7773 A emission line, which was completely 
missing in the old version of the code (e.g. 'Mazz ali et al.l 
2007b). The improved treatment of ionization was necessary 
in order to provide a more accurate estimate of the electron 
density so that recombination rates can be reliably com- 
puted. While ionization rates were obtained from a simple 
analytical form in the older version of the code, in the new 
version they are calculated taking into account the atomic 
and electronic loss processes o f non-thermal electrons explic- 
itly (|Maurer fc Mazzalil [20T0I '). 

Using the i mprov ed treatment of ionization 
(|Maurer fc Mazz"aiil l2010t ). oxygen and other light ele- 
ments are more ionized than they are using the old, simpler 
treatment, while the ionization rates of Fe-group elements 
do not change significantly. Therefore the ionization balance 
of SNe la (which are dominated by Fe-group elements) is not 
infiuenced much by the new treatment of ionization, while 
the ionization balance of SNelb/c, which are dominated by 
carbon and oxygen, is. 

The increased ionization of oxygen and other light el- 
ements results in a larger electron density, which in turn 
favours recombinati on and the formation of th e O I 7773 A 
line, as explained in iMaurer fc Mazzalil (|2010l '). In practice. 
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this turns out to reduce the need for clumping with respect 
to the old version of the code. Using the old ionization treat- 
ment it was possible to mimic the effect of higher ionization 
(higher electron densities) introducing a clumping factor. 
This can explain why clumping seemed to be necessary in 
SNelb/c but not in SNela. As an example, the spectra of 
SN 2007gr can be reproduced using C, — 2.5. This is sufficient 
to guarantee that Fe ill recombines to Fe ll (in SNelb/c 
[Fe III] lines are not seen, unlike SNela). In fact, with the 
new code the Fe lll/Fe ll ratio is similar to the value obtained 
with the old code. Although clumping is not necessarily re- 
quired to obtain good fits to the forbidden line observations, 
it may nevertheless be important e.g. for th e formation of 
the O I 7773 A line l|Maurer fc Mazzalillioiol '). 

The abundances of most elements in SN2007gr (and 
other SNelb/c) are determined from neutral or singly 
ionised species (e.g. C I, O I, Mg I, Na I, Si I, S I, Fe ll, Co ll). 
Higher ionization weakens the flux in these lines. This must 
be offset by somewhat increasing the mass of these elements 
in the nebula. Since the less ionized species remain the dom- 
inant ones even with the updated treatment of ionization, 
the effect on the mass estimate is small. One exception is 
calcium, for which the abundance is obtained from the more 
highly ionised species, Ca ll. In this case, higher ionization 
results in a smaller calcium mass (by a factor of 2-4 depend- 
ing on the individual case). 

The results for day 172 and 390 are shown in Figures 
1, 2, and 3 as red lines. While the fits are similar to those 
obtained with the old model, the masses of some elements 
change somewhat. We now discuss the two epochs in turn. 



4.2.1 Day 172 

The new one-zone model for day 172 is characterised by the 
same outer velocity (5000km s"^) and ^'^Ni mass (0.076Mq) 
as the old code, and the results are similar. However, as we 
discussed above, less clumping is required ((" — 2.5) since 
the ionization degree is higher. In the old model the average 
ionization degree was ~ 1/30 (O ll/O I, Fe lll/Fe ll, etc.), 
the electron density Ue ~ 2.5 ■ 10^ cm"'' and the electron 
temperature Te ~ 5100 K. With the improved treatment, 
the ionization of certain elements (O, Mg, Si, S) is higher, 
^ 20%. The resulting increase in ris favours recombination of 
Fe III, Co III and Ni ill, reducing the need for clumping. The 
ionization of Ca increases less than for lighter elements, and 
it is still only ~ 5%, while that of Fe is essentially unchanged. 
The electron density is now ne ~ 7 x 10^ cm"'' and the 
electron temperature Te ~ 4900 K. 

The increased fraction of highly ionized material weak- 
ens the flux in lines of ions such as O I, Mg I, Si I, making 
it necessary to include more mass in the ejecta. The total 
mass is now ~ 1.8A/q within 5000 km s^^. Oxygen is still 
the most abundant element, with a mass of ~ I.AMq, fol- 
lowed by C (O.25M0). The masses of the other elements are 
comparable to the old model, since their degree of ionization 
does not change, except for Ca, which is reduced by about 
a factor of 3 (O.OI6M0), in accordance with the increase in 
ionization. Calcium is the element for which the new treat- 
ment of ionization has the biggest impact, as the Ca mass 
is estimated from the less abundant singly ionised species. 
Iron-group lines are also mostly from singly ionised species 



(Fe II, Co 11), but these are actually the highly abundant 
lower ionization states. 

The improvements to the synthetic spectrum intro- 
duced by the updated model are subtle, but one can see 
that the recombination line O I 7773 A is now present (Fig- 
ure 1, red line). The mismatch with respect to the data may 
be due to background subtraction: the strength of the line 
seems to be reproduced reasonably well. Alternatively, sev- 
eral Fe II lines near 7800 A may be poorly reproduced if 
their collision strengths are not accurate. With the higher 
ionization, a weak [S 11] line is also visible near 4000 A. 

The sharpness of the peaks of most emission lines con- 
trasts with the more rounded shapes of the synthetic pro- 
files, which are assumed to be parabolae as expected from a 
homogeneous, homologously expanding sphere. This can be 
improved upon using a multi-zone approach. 

4.2.2 Day 390 

As for the day 172 model, the outer velocity (4000km s^^) 
and ^®Ni mass (0.076Mq) used were the same with the new 
code as with the old one. The same clumping as on day 172 
((^ — 2.5) was used. In the optical, the main improvement 
is that now O I 7773 A is reproduced (Figure 2, red line). 
Again because of the higher average ionization, a somewhat 
larger mass is required. The mass within 4000 km s~^ is now 
~ 1.3Mq. The oxygen mass is ~ I.OAIq and the carbon 
mass is 0.14Mq. All other elements have masses similar to 
those used for the model with the old code, except for Ca, 
which is reduced to ~ O.O13M0. The model now has Tie ~ 
4 X lO'^ cm"^ and Te ~ 3400 K. 

In the IR, the quality of the fits with the two versions of 
the code is very similar (Figure 3, red line). As with the old 
code, a rather small mass of silicon is sufficient to reproduce 
the [Si l] line near 1.6/x. 



5 SHELL MODELS 

In this section we present shell models for nebular obser- 
vations of SN 2007gr at 172 and 390 days after explosion. 
The models were computed using the updated version of the 
code. The results of one-zone models (in particular abun- 
dances and their ratios) were used as the starting point 
for the shell models, which were however optimised inde- 
pendently. In contrast to the one-zone modelling we do not 
present 'best fit' models for each epoch, but we try to find 
one model which can reproduce the observations at 172 and 
390 days simultaneously. Only sodium and magnesium mass 
estimates change between the two epochs. 

While there is no IR spectrum at 172 days, there is an 
IR observation of SN 2007gr at 190 days after explosion. 
In order to be able to merge the optical and the IR into 
a single spectrum, and lacking any IR photometry on day 
172 (Hunter ot al 2009), we extrapolated the 190-day IR 
observation to 172 days. To do this, we computed the IR-to- 
optical flux ratio of synthetic models at 172 and 190 days, 
and scaled the 190-day IR observations to reproduce the ex- 
pected ratio at 172 days. Since the difference is 18 days only, 
this extrapolation should be reasonable. We can then com- 
pare the optical and IR flux at both 172 and 390 days. Shell 
models for the day 172 optical spectrum and the rescaled IR 
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Figure 4. SN 2007gr: models for the optical spectrum at day 172 
(29 Jan 2008, black) with the shell version of the nebular code, 
which includes an improved treatment of ionization and oxygen 
recombination (green). 

spectrum of day 190 are shown in Figures 4 and 5, respec- 
tively, while fits for the day 390 optical and IR spectra are 
shown in Figures 6 and 7, respectively, as green lines. 

The total mass of our model is ~ l.OM© below 
6000 km s~^ (the model becomes inaccurate at higher veloc- 
ities because 7-ray deposition is very small at low density). 
This is in rough agreement with the one-zone models (see 
Tables 1 and 2). In general, the masses derived with the shell 
model are slightly lower. This is expected since one-zone 
models usually produce too much flux around the peak of 
the lines, which leads to an overestimate of the line luminos- 
ity and thus of the amount of emitting material. Mass ratios 
of individual elements are similar to the one-zone models. 

One drawback of the nebular modelling approach is the 
poorly constrained distribution of ^®Ni. The one-zone model 
assumes that ^''Ni is distributed homogeneously. Some con- 
straints on the precise distribution of ^''Ni can be derived 
from the profile of the Fe 11 blend near 5200 A, but this 
is usually weak in SNelb/c. In addition, not all ^®Ni dis- 
tributions can consistently reproduce the evolution of the 
spectrum from 172 to 390 days. In our shell model ^®Ni is 
more centrally concentrated than in the one- zone models, in 
agreement with theoretical expectations. The uncertainty on 
the total mass caused by this assumption should be ~ 10%. 

We use a clumping factor ^ = 2.5, which gives the best 
overall agreement with the observations, although individual 
lines or line ratios can sometimes be better reproduced using 
different values. Any moderate value, ^ < 10 seems possible, 
taking into account uncertainties in the atomic data and the 
calculation procedure. 

Since the same model is used for both epochs (apart 
from differences in the Na and Mg content), we discuss the 
two epochs together, focussing on individual elements. 
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Figure 5. SN2007gr: models for the IR spectrum at day 190, 
rescaled to day 172 (16 Feb 2008, black) with the shell version of 
the nebular code, which includes an improved treatment of ioniza- 
tion and oxygen recombination (green). Lines marked in red are 
tentative identification of features that are not reproduced in the 
current version of the code: they could be mostly recombination 
lines of Mg I and O I. 



The carbon distribution is mainly derived from mod- 
elling the [C l] 8727 A emission at 172 days. The flux in this 
region is a blend of [C l] and Ca ll emission and the ob- 
served flux profile cannot be reproduced by any model with 
our nebular code. It may be that the Ca ll fiux is scattered 
by the carbon line, a process which is not simulated by the 
code. Alternatively, [Fe ll] lines may be stronger than pre- 
dicted by the model, as discussed in Section 4.2.1. At day 
172, we try to reproduce the fiux at the wavelength of the 
[C l] line. Using the same carbon abundance, at 390 days 
the synthetic [C l] 8727 A line is stronger than in the obser- 
vations. This may be an indication that at day 172 the line 
was powered by scattering of Ca ll fiux: since this process 
is not treated in our code, we may have had to assume a 
higher carbon mass than is really present in the SN ejecta. 
Since no other line of carbon is visible, the estimate of the 
carbon mass may thus be uncertain by a factor of ~ 2. 

The oxygen distribution and mass are mainly derived 
from modelling the [O l] 6300, 6363 A doublet. The observa- 
tions are matched rather well by the synthetic flux at both 
epochs. The strength of the [O l] 5577 A line depends on 
clumping. Using a clumping factor = 2.5 it is reproduced 
well at 172 days. At 390 days the line is no longer visible, 
since its formation is strongly favoured at higher densities. 
With the new code, several permitted O I lines which arise 
from recombination can be modelled. The O I 7773 A line is 
not reproduced with sufficient strength at 172 days, but at 
390 days the synt hetic ffux is in goo d agre ement with the 
observations (see iMaurer fc Mazzalil ||2010D for a detailed 
discussion of the formation of this line) . The O I 8447 A line 
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Table 2. Results of Shell Models. The masses listed here are located below 6000 km s^^. About 0.4 Mq 
of oxygen are expected above 6000 km s^^. Since there is no information about the helium content of SN 
2007gr the values for total mass are lower limits. 



epoch 


Mass 


M(C) 
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M(Na) 


M(Mg) 


M(Si) 


M(S) 


M(Ca) 


M{56Ni) 


days 


M0 


Mq 


Mq 


Mq 


Mq 


M© 


Mq 


Mq 




172 


1.02 


0.15 


0.76 


5.0e-5 


9.5e-4 


0.005 


0.020 


0.009 


0.076 


390 


1.05 


0.15 


0.76 


5.0e-4 


0.028 


0.005 


0.020 


0.009 


0.076 
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Figure 6. SN 2007gr: shell model for the optical spectrum at day 
390 (3 Sept 2008, black) with the shell version of the nebular 
code, which includes an improved treatment of ionization and 
oxygen recombination. The green line shows a model computed 
for clumping = 2.5, the blue line is a model computed for C, = 66. 



Figure 7. SN 2007gr: shell model for the IR spectrum at day 
390 (3 Sept 2008, black) with the shell version of the nebular 
code, which includes an improved treatment of ionization and 
oxygen recombination . The green line shows a model computed 
for clumping = 2.5, the blue line is a model computed for = 66. 



is blended with strong Ca ll emission and cannot be identi- 
fied. The O I 9264 A line is too weak at day 172. Since its 
formation is directly related to the 7773 A line this is ex- 
pected. In the IR, the O I 1.129 /xm line is swamped by [S l] 
1.131 ^im. At 390 days the synthetic flux of the O I 7773 A 
and 9264 A lines reproduces the observations, taking into 
account the possible background. The O I 8447 A line still 
cannot be identified in the observations. 

Sodium and magnesium are modelled fitting the Na I D 
and Mg l] 4571 A observations, but different masses are re- 
quired at different epochs (see Sec. 4.2), which indicates that 
these ions are not accurately treated. In fact, the IR spec- 
trum at day 190 shows a number of emission lines which are 
not reproduced by the code (see red marks in Fi gure 5). The 
lineat 1.5 /i was identified as Mg i 1.503^ by IPozzo et al.l 
Several other lines can be identified as permitted 
Mg I transitions, as indicated in Fig. 5. These lines could be 
produced by recombination of Mg ll. To our knowledge, this 
is the first time these lines are identified in a nebular spec- 
trum of a SN. The fiux in these lines is intrinsically weak 
(note the factor ~ 10 difference in the vertical scale of Fig. 5 



with respect to Fig. 4). Inclusion of a detailed treatment of 
Mg 11 recombination may result in a better fit of these lines 
and in a more accurate determination of the Mg mass. 

Two other lines are not identified in the IR spectrum at 
day 190. The line at 1.3 /.im may be O I 1.315 /im. This would 
be part of the O ll recombination process. The emission 
near 1 .43 /im was identified as [Fe l] 1.443 /^m bv lPozzo et al.l 
(l2006l ) in the spectrum of a SN la obtained very early in the 
nebular phase. While this identification may be possible, 
the lack of any other [Fe l] lines, in particular 1.294 and 
1.355 fj,m, which arise from the same lower level (the ground 
state of Fe l) and should have similar strength is a strong 
argument against this option, which remains however the 
only one we could find. 

The silicon distribution is modelled fitting the [Si l] 
1.588, 1.608 and 1.646 /xm lines at 390 days, when the [Si l] 
feature can be observed clearly. The line is blended with 
[Fe ll] lines, and the same procedure described in Section 
4.1.2 is used to determine the mass of silicon. At 172 days 
this feature is blended with [Co ll] emission. Also with the 
shell model the silicon mass remains remarkably small. 
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The sulphur distribution can be obtained modelling the 
[S ll] 4070 line at 172 days, when it is more clearly visible, 
as well as the [S l] 1.083 and 1.131 ^m lines. The [S l] 7722 A 
line is blended with O I 7773 A. Since it is difficult to iden- 
tify individual lines in the region around 4000 A, we fix the 
sulphur mass by fitting the IR lines accurately, while aim- 
ing at an acceptable agreement for the 4070 A line. At 390 
days the synthetic IR lines seem to be too weak, but the 
observations are very noisy in that spectral region. 

The calcium distribution and total mass are rather un- 
certain. The ratio of Ca ll] 7300 A and the IR-triplet de- 
pends sensitively on density and clumping. In addition, the 
IR-triplet may interact with C I 8727 A, a process which is 
not simulated in the code. We fix the Ca mass modelling the 
7300 A emission at 172 days. This gives a consistent result 
also at day 390. 

Since we cannot obtain an estimate o f the helium 
mass, the total mass in SN 2007gr is uncertain. iHunter et al.l 
l|2009h derive a total mass of ~ 2 Mq from light curve mod- 
elling, which would imply ~ 0.5 M0 of He. 

Finally, the spectra for day 390 show an apparent mis- 
match between the optical and IR Fe lines: most IR lines 
seem to be overestimated while the [Fe ll] emission near 
5200 A is too weak. In order to test whether this may be the 
effect of clumping, we have computed a model with a clump- 
ing factor which brings the optical synthetic spectrum into 
reasonable agreement with the observations. This requires 
extreme clumping ((^ — 66): the model is shown in Figures 6 
and 7 as a blue line. Increased clumping leads also to a re- 
duction in the flux of the [C l] lines in the optical. However, 
despite the large change in the fiux of the optical [Fe ll] lines, 
the IR [Fe ll] lines are almost unaffected: increased clumping 
leads to a reduction of the population of the highly excited 
levels from which the optical lines are formed, but leaves 
the much larger population of the less excited levels from 
which the IR lines form basically unchanged. Another pos- 
sibility would be to increase the ^®Ni mass. The [Fe ll] com- 
plex near 5200 A can be reproduced in sufficient strength 
for M(^'^Ni)= O.IM0, but then the IR Fe fines become too 
strong. Further increase of the clumping seems unlikely: al- 
ready with ( — 66 the spectrum at day 172 is poorly re- 
produced. Considering the uncertainties in the background 
subtraction, the Fe atomic data, the ^®Ni distribution and 
the geometry, it is not possible to determine the degree of 
clumping accurately, at least using the code in its present 
state. 



6 DISCUSSION 

The availability of a simultaneous optical-IR spectrum for 
SN 2007gr makes it possible to attempt an accurate esti- 
mate of the ejecta mass via a census of all elements that are 
expected to be abundant in the ejecta of a SNIc. 

Our models indicate that SN2007gr ejected ~ IMq of 
heavy elements. Our models are sensitive to material ex- 
panding with V < 6000 km s~^. We adopt here the values 
from the shell model, since this offers the best reproduction 
of the observed line fluxes and proflles. We use a ^^Ni mass of 
0.076 Mq, as determined f rom the luminosity at light curve 
peak (|Valenti et al.ll2008l : [Hunter et al.ll2009l '). and confirm 
that that value is a reasonable estimate. We find that oxygen 



is the dominant constitutent of the heavy-element ejecta, 
with a mass of ~ O.8M0. The IME abundances, many of 
which can be determined thanks to the IR spectrum, are 
quite small. In particular the Si mass, which is determined 
via the [Si l] 1.6 /im line, is only O.OO5Af0. Since the Si mass 
is estimated from neutral silicon, which is the dominant ion, 
this result is not much infiuenced by ionization. 

Our estimate of the Ca mass carries a larger uncertainty 
than for other elements since it is estimated from lines of 
Ca 11, which is less abundant than Ca I. Uncertainties in 
the ionization treatment can affect the mass estimate signif- 
icantly (a factor of 3 or so). 

Although we use moderate clumping — 2.5), it is not 
possible, at the current level of accuracy, to determine the 
exact degree of clumping. 

The rather low value of the ejected mass may be in- 
creased if the mass of He and possibly Ne, which may 
be present t o som e extent, could also be accounted for. 
iHunter et al.l (|2009l ) estimated an ejected mass of ~ 2 M© 
from light curve models. If this estimate is correct, the He 
mass of SN2007gr should be ~ 0.5 - l.OM©. This may be 
consistent with early-time observations of He I absorption, 
depending on the degree of mixing of He and radioactive 
material in the outer layers of the SN ejecta. 

The result that the mass of the CO core ejected by 
SN 2007gr was ~ IMq suggests that the SN progenitor was 
a star of comparatively low mass, possibly similar to the 
progenitor of the prototypical low-energy SNIc 19941. In 
the case of SN 1994 1 a star of ~ ISM^ was proposed a s 
the likely progenitor (|Nomoto et al.lll994l : [Sauer et al.ll2006l '). 
The star would have lost its outer hydrogen and helium lay- 
ers through binary interaction. The progenitor of SN 2007gr 
may have had a similar evolutionary history, with less strip- 
ping of the outer layers. The estimate of the progenitor mass 
is also bro adly consistent with t he observational constraints 
derived by Crockett et alj l|2008l ) from pre-explosion images. 

Hunter ct al. (2009, ') estimated a '^'^Ni mass of 0.076 Mq. 
Using this value, we can reproduce the temporal evolution of 
the nebular observations from 172 to 390 days. At 172 days 
the strength of the [Fe 11] blend near 5200 A is reproduced 
reasonably well, while at 390 days the synthetic flux at these 
wavelengths is too low. Uncertainties in a number of para- 
menters, like the degree of clumping, the ^^Ni distribution, 
Fe collision strengths, geometry, reddening and background 
subtraction could be responsible for this. A high clumping 
factor ~ 60) or a larger ^'^Ni mass (~ O.IMq) can approx- 
imately reproduce the strength of the [Fe 11] blend. However, 
an increased ^^Ni mass leads to an over-estimate of less en- 
ergetic [Fe 11] lines redwards of 6000 A. Also, a ^^Ni mass of 
more than 0.09Mq is inconsistent with the 172 days spec- 
trum. Therefore we estimate a total ^®Ni mass between 0.07 
and O.O9M0, but favour the lower values. Our results are in 
fact consiste nt with the estimate s of the total and ^^Ni mass 
presented in IHunter et al.l l|2009l ). 

The elemental abundances that we derived are gener- 
ally consistent with the nucleosyn thetic yields of a « ISMq 
star at roughly solar metallicity (IWoos lev fc Weaver! Il995l : 
[Thielemann et al ] |l996l : lNomoto et al.l ll997. 200^). The ma- 
jor discrepancy with respect to those results is actually for 
silicon, which we flnd to be about one order of magnitude 
less abundant. Our results are unlikely to be grossly in er- 
ror because of an incorrect estimate of the ionization: Si I, 
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the ion which is observed in the IR and from which the Si 
mass is determined, is the dominant silicon ion. On the other 
hand, possible uncertainties on the collision strengths of Si I, 
or of the Fe ll lines which also contribute to the feature at 
1.65 fim from which the silicon abundance is measured, may 
cause uncertainties. If the [Fe ll] lines were completely re- 
moved, and the collision strength for the [Si l] lines remained 
the same, the silicon mass would increase by a factor of ~ 3. 
At the same time, the results cited above, while in over- 
all agreement, differ the most for IME (Si, S, Ca), and more 
recent calculation (S.E. Woosley, priv. comm.) show a some- 
what smaller silicon mass, although again a factor of 10 error 
is unlikely. The most likely reason for the discrepancy may 
lie with the collision strengths. A new calculation of these 
values would be desireable. 

The detection of radio emission from SN2007gr, indi- 
cating the presence of material moving at s ub-relativistic 
velocities (i; ~ 0.6 c) in two opposite jets (jParaei et al.l 
[2010) is a very interesting finding. Normally, SNe associ- 
ated with relativistic outflows such as GRBs or XRFs have 
exceptionally large kinetic energies and massive progenitors. 
This does not seem to be the case for SN 2007gr. Actually, 
the nebular spectra of SN 2007gr show that the core of the 
SN has rathe r low velociti e s, in agreement with the early 
phase, where IValenti et all (|2008h found ma^ximum veloci- 
ties of « 10000 km s~^. Actu ally, in a paper that ap peared 
after submission of this paper, ISoderberg eT^ (|2010l ) argue 
that the radio data are consistent with a non-relativistic ex- 
plosion, in line with our optical results. 

Another peculiarity of SN 2007g r was the detectio n 
of strong carbon lines at early times (|Valenti et al.ll2008l ). 
These lines were observed with velocities between 7000 and 
10000 km s~\ In the nebular spectra the [C l] 8727 A line 
is detected, suggesting the presence of carbon also at low 
velocities (below 5000km s^^). Modelling indicates a C/0 
ratio of ~ 0.1 — 0.2, which is consistent with the inner 
layer s of the CO core of a massive star ([Wooslev et al.l 
19941: [Thielemann et al] 1 19961 : iNomoto et all Il997l . l2006l : 



Maeder fc Mevnetll2003l . l2005l ) 



Finally, GRB/ SNe are thought to be very aspherical 
iMaeda et aLllioOq ). Some evidence that all SNcIb/c are as - 
pherical, to some exten t, was found bylMaeda ct al] (|2008l ) 
using late-time spectra. iMaurer et al.l (|2010l ) confirmed this, 
but suggested that for SNelb/c which are not particularly 
energetic there may be only a rather weak correlation be- 
tween the velocities of the inner ejecta ("the core") as de- 
rived from nebular spectroscopy and those of the outer 
ejecta, as inferred from the spectra near maximum. The 
highest velocities detected in the optical spectra {v ~ 0.1c) 
are always much lower than the sub-relativistic velocities 
which can give rise to radio emission, so SN 2007gr may be 
a case where a very small amount of ma terial was ejected a t 
V ~ 0.6 c in an otherwise normal SNIc JParagi et al.l 201ol). 
Howev er, the recently published results of ISoderberg et ed] 
l|2010l ). who rule out material moving at w > 0.2c are in 
much better agreement with our results. From the nebu- 
lar spectra we find no evidence of major core asphericities 
in SN 2007gr. This does not rule out that so me asphericity 
may affect the innermost regions of the ejecta (iMaurer et al.l 
I2OIQ V This may also help with the missing silicon if more 
fallback occurs in the less nuclearly processed equatorial re- 
gions. Clearly, further study of these very interesting events 



is necessary before a coherent picture that encompasses the 
entire range of SNelb/c properties and their relation to the 
those of the progenitor star can emerge. 
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